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Molecular analysis of an odorant-binding protein gene 
in two sympatric species of Lutzomyia longipalpis s.l. 
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Lutzomyia longipalpis s.l. is the main vector of American visceral leishmaniasis (AVL) and occurs as a species 
complex. DMA samples from two Brazilian sympatric species that differ in pheromone and courtship song produc- 
tion were used to analyse molecular polymorphisms in an odorant-binding protein (obp29) gene. OBPs are proteins 
related to olfaction and are involved in activities fundamental to survival, such as foraging, mating and choice of 
oviposition site. In this study, the marker obp29 was found to be highly polymorphic in Lu. longipalpis s.l., with no 
fixed differences observed between the two species. A pairwise fixation index test indicated a moderate level of ge- 
netic differentiation between the samples analysed. 
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The sandfly Lutzomyia longipalpis s.l. is the main 
vector of American visceral leishmaniasis (AVL) and oc- 
curs as a complex of cryptic species (Bauzer et al. 2007, 
Araki et al. 2009). A number of integrative studies using 
different approaches, such as molecular polymorphisms 
(Bauzer et al. 2002a, b, Bottecchia et al. 2004, Araki et 
al. 2009, Lins et al. 2012), microsatellites (Maingon et 
al. 2003), pheromones (Hamilton et al. 1999a, b, Souza 
et al. 2004) and courtship song analysis, have revealed 
the existence of at least five species in Brazil (Souza et 
al. 2002, Araki et al. 2009). Lu. longipalpis s.l. males 
exhibit morphological polymorphism with regard to the 
number of pale abdominal tergal spots and can present 
just one pair of pale spots on the fourth abdominal terg- 
ite (IS phenotype) or two pairs of pale spots, one on the 
third and the other on the fourth abdominal tergite (2S 
phenotype). Intermediate phenotypes (a small spot on 
the 3rd tergite in addition to the spot on the 4th tergite) 
are observed in high frequencies in some localities, indi- 
cating an intraspecific polymorphism (Ward et al. 1988). 
In contrast, intermediate forms are rare or nonexistent 
in localities where two Lu. longipalpis s.l. cryptic spe- 
cies occur in sympatry (Bauzer et al. 2002a, Araki et 
al. 2009). These tergal spots contain pheromone glands 
that are associated with Lu. longipalpis sexual commu- 
nication (Lane et al. 1985). In the Brazilian locality of 
Sobral, state of Ceara, a sympatric species identified by 
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the 2S phenotype was shown to produce the cembrene-1 
(C, 0 )-type pheromone, whereas a species associated with 
the IS phenotype was shown to produce the 9-meth- 
ylgermacrene-B (C ]6 )-type pheromone (Lane et al. 1985, 
Hamilton et al. 1999a, b). 

Proteins related to olfaction are involved in activi- 
ties that are fundamental to survival, such as foraging, 
courtship, mating and choice of oviposition sites (Hallem 
& Carlson 2004). Among these molecules, odorant-bind- 
ing proteins (OBPs) have been described as important 
components in the recognition of odours (Hekmat-Scafe 
et al. 2002) and several OBPs from different insect spe- 
cies have been cloned and sequenced (Xu et al. 2003, 
Zhou et al. 2004, 2008). Furthermore, an increasing 
number of OBPs have been identified in several non- 
sensorial tissues, such as salivary glands (Abdeladhim 
et al. 2012), head and body (Li et al. 2005, Gonzalez-Ca- 
ballero et al. 2013), reproductive organs (Azevedo et al. 
2012), fat bodies and seminal fluid (Liu et al. 2010, Sirot 
et al. 2011). The wide distribution of these proteins in the 
organism suggests that they can perform several other 
physiological functions (Pelosi et al. 2006, Pelletier & 
Leal 2009). Phylogenetic analyses and the distribution 
and orientation of OBP genes in insect genomes have 
revealed evidence of a complex series of duplication and 
rearrangement events, facts that suggest that this gene 
family evolves rapidly (Hekmat-Scafe et al. 2002, Vieira 
et al. 2007) and are therefore potentially good markers 
for use in population genetic studies. 

In this study, we analysed polymorphisms in the obp29 
gene to characterise the molecular variation of sympatric 
species of Lu. longipalpis s.l. and to hypothesise a possi- 
ble role of the encoded OBPs in the reproductive isolation 
and adaptation of this important insect vector. 

The genomic DNA from samples collected in the 
Brazilian locality of Sobral (3°41'S 40°20'W) was the 
same as that used by Bauzer et al. (2002a.) Genes that 
were isolated from Lu. longipalpis s.l. and searched 
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for using the Antennae EST database developed in our 
laboratory (Dias 2008) directed the choice of specific 
oligonucleotide primers for the obp29 gene. Polymer- 
ase chain reaction (PCR) reactions were performed ac- 
cording to the Go Taq DNA Polymerase instructions 
manual (Promega) using a thermocycler (GeneAmp 
PCR System 9700/Applied Biosystems) with forward 
(5'-GGGAGTGAGAAAAATGAGATCAAA-3') and 
reverse (5'-ACTTGACATTGCTTTTCTGTGCAGG-3') 
OBP29 primers. The following programme was used: 30 
cycles of 95°C for 30 sec, 60°C for 30 sec and 72°C for 
1 min. The obtained fragments were purified using the 
QIAquick PCR Purification kit (QIAGEN) and cloned 
using the pGEM-T Vector kit (Promega). Plasmid DNA 
was isolated using the mini-preparation alkaline lysis 
method (Sambrook & Russel 2001) in 96-well MicroW- 
ell Plates and the DNA was cleaned using filter plates 
(Multiscreen-Millipore). Sequencing of the cloned frag- 
ments was carried out using the ABI Prism Big Dye Ter- 
minator Cycle Sequencing Ready Reaction v.3.1 kit (Ap- 
plied Biosystems) and the ABI 3730 Sequencer by the 
Program on Technological Development of Health (Fi- 
ocruz, Rio de Janeiro, Brazil) sequencing service (Otto 
et al. 2008). The obtained sequences were edited using 
the Bioedit v.7.0.5 programme (Hall 1999) and an align- 
ment was performed using the CLUSTALW full multiple 



alignment option within the BioEdit software. Different 
measures of genetic polymorphism were calculated us- 
ing DnaSP v.5 (Librado & Rozas 2009), including the 
number of segregating sites (S), the nucleotide diversity 
pi (?i) which is the average number of nucleotide sub- 
stitutions per site between any two sequences and the 
parameter theta (8), which is estimated from the number 
of segregating sites. The estimated n and 9 parameters 
were used to perform Tajima's D test of neutrality (Taji- 
ma 1989). The Hudson, Kreitman & Aguade neutrality 
test (HKA) was performed using the period gene (Bau- 
zer et al. 2002a) as the second locus. Both neutrality tests 
were executed using the DnaSP v.5 software (Librado & 
Rozas 2009). A population structure analysis was con- 
ducted by separating two groups of sequences based on 
the abdominal spot phenotype (IS or 2S). The population 
subdivision was calculated using the programme ProS- 
eq3.0 (Filatov 2002) and the significance of the pairwise 
fixation index test was estimated by a permutation ap- 
proach. The sequences were submitted to GenBank (ac- 
cessions KF669570-KF669605). 

A 685-bp opb29 fragment encompassing the coding 
region and part of the 3' UTR was obtained. A total of 
18 consensus sequences were analysed from the Sobral 
IS and 2S populations. Of the 685 analysed nucleotides 
sites, 93 (13.6%) were variable and the nucleotide diversi- 
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Alignment of the variable sites in the Lutzomyia longipalpis s.l. odorant-binding protein 29 haplotypic sequences. Dots indicate nucleotides shared 
with the first sequence. Positions of the polymorphic sites are indicated on the top of the figure. Non synonymous sites are shown in gray. 
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ty estimation was similar for the two species analysed (ti 
= 0.035 for Sobral IS and n = 0.0331 for Sobral 2S). Simi- 
larly, the parameter 8 was found to be essentially the same 
for the two species (0.0310 and 0.0309 for Sobral IS and 
Sobral 2S, respectively). Figure shows an alignment of 
all the polymorphic sites observed in the analysed opb29 
fragment. Twenty three polymorphic sites were exclusive 
to Sobral IS, 24 were exclusive to Sobral 2S and 54 were 
shared between the species; non-synonymous mutations 
were found at 12 sites (6 exclusive to SIS, 3 exclusive to 
S2S and 3 shared). A moderate and significant genetic 
differentiation (Fst = 0.1098; p < 0.001; 1,000 permuta- 
tions) was computed between the two species. The genet- 
ic divergence and polymorphism data were used to test 
departures from neutrality, with both Tajima's D statis- 
tics (Tajima 1989) and the HKA test (Hudson et al. 1987) 
indicating no departure from neutrality. 

The performed molecular analysis revealed no fixed 
differences between the two analysed populations. Al- 
though some exclusive synonymous and non-synony- 
mous mutations were found, a clearer pattern of differ- 
entiation was not obtained. Therefore, we cannot infer a 
possible role of the obp29 gene in responses to specific 
pheromones or in the process of reproductive isolation. 
The level of genetic differentiation was not as high as 
that observed for other nuclear markers (Bauzer et al. 
2002a, b, Bottecchia et al. 2004, Araki et al. 2009, Lins 
et al. 2012). This finding can be explained by the recent 
origin of these species within the Lu. longipalpis com- 
plex and the consequent retention of ancestral polymor- 
phisms. Alternatively, one might also consider introgres- 
sion events if reproductive barriers allowed rare events 
of hybridisation. Future studies expanding this analysis 
to other obp genes might help in the elucidation of a pu- 
tative functional role of OBPs in the speciation process 
occurring mLu. longipalpis s.l. 
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